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ABSTRACT 

O 
O 

CN ■ The Letter takes up a question of what radio emission is produced by electrons at 

the very acceleration site of a solar flare. Specifically, we calculate incoherent radio 
emission produced within two competing acceleration models — stochastic accelera- 
tion by cascading MHD turbulence and regular acceleration in collapsing magnetic 
traps. Our analysis clearly demonstrates that the radio emission from the accelera- 
|Vj tion sites: (i) has sufficiently strong intensity to be observed by currently available 

radio instruments and (ii) has spectra and light curves which are distinctly different 
■ in these two competing models, which makes them observationally distinguishable. 

In particular, we suggest that some of the narrowband microwave and decimeter 
continuum bursts may be a signature of the stochastic acceleration in solar flares. 

a 

Subject headings: Sun: flares — acceleration of particles — turbulence — diffusion- 
Sun: magnetic fields — Sun: radio radiation 

> ■ 
ON 

in 

' 1. Introduction 

t> 

Acceleration of charged particles is an internal property of energy release in solar flares, 
which has not yet been fully understood in spite of a significant progress achieved recently (e.g. 



Aschwandenl 120021 : IVilmer &: MacKinnon! 120031 ). A traditional way of getting information on 



^ . the accelerated electrons in flares is the analysis of the hard X-ray (HXR) emission produced 
by nonthermal bremsstrahlung. However, because the bremsstrahlung intensity increases with 
the density of the ambient plasma, it is likely that in most cases the acceleration site and 
HXR emission site are spatially separated; therefore, the HXR emission does not carry direct 
information on the acceleration site. This concept of distinct acceleration, pro pagation, and 
emiss ion regions was then inherited by solar radio astronomy (e.g., Fig. 9 in iBastian et al. 



19981 ). which looks relevant to relatively weak events of electron acceleration visualized by 
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coherent emission of type III groups and of accompanying metric spikes (e.g., Fig. 10 in 
Bastian et a"Dll998h . 



However, it is well known that a charged particle produces electromagnetic emission as 
it moves with acceleration. Stated another way, fast electrons must produce radiation imme- 
diately at the acceleration region with intensity and other characteristics defined by type of 
the acceleration, or more precisely, by the type of fast electron trajectories in the acceleration 
region. We show in this Letter that typically this emission has a spectral peak at the microwave 
range, which makes the radio observation the most suitable to study the acceleration region in 
flares. 

By now, a huge variety of acceleration mechanisms and models has been proposed and 
develope d for the solar flares. Acceleration by DC electric fields, both sub-Dreicer and super- 
Dreicer, jHolmanl Il985l : iTsunetal Il985l : lHolman^_Benkalll992l; iLitvinenkol Il996l. boool. I2OO3I) : 



Mil ler et al.lll996l : lMillerill997l : IPark et al.lll997k IPrvadko & Petrosian 



stochastic acceleration by turbulent waves ( Hamilton fc Petrosian 1 19921; 



Ferm i ) acceleration in collapsing magnetic traps (Somoy fc Kosugil 119971 ; ISomov fc Bogachev 



20031 : iKarlickv fc Kosugil booi iBogachev fc Somovi [20051 . 120071 . I2OO9I ): all are currently consid- 



199 



Petrosian et al 



1994 



Bykov fc Fleishman 



20091 ). the classical diffusive shock acceleration ( Aschwanden 2002 ); the regular (betatron plus 



ered in the context of solar flares. 

To illustrate potential ability of radio observations to detect the radiation from the flare ac- 
celeration site and to distinguish then between competing acceleration mechanisms, we calculate 
here radio emission generated within two distinct acceleration models — stochastic acceleration 
by a turbulence spectrum and regular acceleration in collapsing traps. Radio emission of flares 
is known to be produced by a variety of emission mechanisms including gyrosynchrotron (GS) 



Bastian et al. 


1998; 


rt! „„„„„ 

Nindos et al. 


2008) 



inate from acceleration sites, while others — from electrons trapped in closed magnetic loops 
or from electrons propagating along open field lines. Based on our analysis, we suggest that 
some of the narrowband microwave and decimeter continuum bursts may be a signature of the 
stochastic acceleration in solar flares, while the collapsing trap acceleration must reveal itself 
in drifting to higher frequencies microwave GS bursts. 



2. Radio Emission from a Region of Stochastic Acceleration 



Basically, various models of stochastic acceleration differ from each other by the accelerat- 
ing agent (the plasma or MHD eigen-mode responsible for the wave-particle energy exchange) 
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and presence or absence of some pre-acceleration (injection) process. To be specific, we assume 
a 'pure' stoch astic acceleration process when the electrons are accelerated directly from the 
thermal pool ( Hamilton fc Petrosian 1992 ; Miller et alill996l ). perhaps as a result of MHD tur- 
bulence cascading towards the smallest scales involved into resonant interaction of the waves 
with thermal or weakly superthermal electrons. 

Within this model the MHD turbulence is created at some large scale and then a broad 
spectrum of the turbulence develops due to the turbulence cascading. As soon as small-scale 
waves capable of resonant interaction with electrons from Maxwellian tail are produced they 
start to accelerate those electrons. This process can be modeled by growing a power-law tail 



(sf., e.g. spectra of accelerated e lectrons presented by lHamilton &: Petrosianlll992l ; iMiller et al. 
19961 : iBvkov fc Fleishman] l2009h 



N(E) = (S(t) - 1) 



w(> E ) • E ( 
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S(t)-1 



exp 



E 



(1) 



E br {t) 

where the time-dependent acceleration is modeled by increasing the break energy E br {t) and 
hardening the energy spectrum (decreasing spectral index S(t)). This nonthermal distribution 
of accelerated electrons is assumed to match the original Maxwellian distribution at a certain 
energy E ; n(> E ) is evidently defined by the matching condition: 



n(> En 
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E 3 
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exp 
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exp 
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where n e and T e are the number density and temperature of the thermal electrons, k is the 
Boltzman constant. Figure [1] shows a sequence of the electron spectra as the electron acceler- 
ation progresses. 

Let us consider the radio emission produced by accelerated electrons with the spectrum 
(CQ) at the acceleration region. We note that gyrosynchrotron (GS) emission by nonrelativistic 
and weakly relativistic electrons available during an initial phase of the acceleration modeled 
by Eq. (CQ) is inefficient; the flux of the GS emission remains ty pically very small un til sufficient 
number of electrons is accelerated up to a few hundred keV[j ( iBastian et al.lll998l ). However, 
along with the regular magnetic field, there is a spectrum of turbulent waves (those accelerating 
the electrons) at the acceleration site. The nonthermal electrons, interacting with those random 
waves, experience spat ial diffusion and so produce so called Diffusive Synchrotron Radiation 
(DSR. lFleishmanll2006l ). which we calculate here. 



1 We note that in case of big fla res, large numbers of GS-producing electrons can already be generated during 
a preflare phase ( Asai et all 120061 ). In such cases we have in mind an even earlier stage of acceleration (e.g., 



Battaglia et al. 20091 ) . when the 100 keV electrons are not yet numerous. 
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Fig. 1. — Distribution function of electrons for different S and the following parameters: the 
thermal electron number density n e = 10 10 cm~ 3 ; the electron temperature T e = 10 6 K; the 
matching energy is Eq = QkT e . The solid curve denotes the original Maxwellian distribution. 



Neglecting for simplicity the plasma gyrotropy we can take the refractive index of the radio 
waves in the form 



Tin 



U) 2 



pe 



Ane 2 n P 



m 



(3) 



Then, spectral and angular distribution of the energy radiated by a relativistic charged particle 
with a given Fourier tran sformed acceleration w,,,/ during time T of the particle motion in an 
external field is given by (ILandau fc Lifshitzlll97ll ) 
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where 



to = to 



(l-^V^R). (5) 
In the nonrelativistic case u/c<l (7 = a/1 — v 2 /c 2 ~ 1) and uo' m oj, Eq. (j3J) reduces to 



£n,u = Ve%- |[n x w^]] 2 , 

c 6 



(6) 
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where Q is the particle charge and n is the unit wave vector of the radiation. Eq. ([6]) shows 

2 

that the radiation in a given direction n is defined by the acceleration componen t |w M i | 



| [n x w^]) 2 transverse to n. Similarly to the derivation in ultrarelativistic case (jFleishman 



20061 ). the transverse component of the acceleration can be expressed via temporal and spatial 



Fourier transform of the external Lorentz force, F° Q 



(27r) 3 f 

I w w± | 2 = / dq dq5(u - q + qv){5 af3 - n a np)F£ A F£ A , (7) 

where M is the mass of emitting particle and V is the source volume. 
For electric component of the Lorenz force F = QE we have 

(Sap - nan^F^F^ = Q 2 (5 aP - n a n p )E^E^ (8) 

where 

TV 

KA= (2^4^(90, q), (9) 
K n n{q (\,c\) is the corr elation tensor of the random electric field, such as J dq dqK aa (qQ, q) = 



(E 2 st ) (|Toptygirj|l985|). 



For magnetic component of the Lorenz force the corresponding expression is different 
(<W - n a n p )F- M F^ = g (v 2 S aP - v a v p - [n x v ]«[n x v]„) ^„ q < q = 

(10) 

Q% (n a n, + ^5 a , - (nv)^+^) B^B^. 
Similarly to Eq. (Q we have 



TV 

W) 

where K a/ 3(q , q) is the correlation tensor of the random magnetic field, such as J dq dqK aa (q , q) 
(B 2 t ). Thus, the DSR intensity, I n uj = S n ^/T, of a nonrelativistic particle in the presence of 
random magnetic field is 

Q 4 v 2 sf . (nv) 2 .v a np + n a vp\ , , 

= V g 2 7rM 2 5 / dq ° dcL °( UJ + ( n » n P H r" ~ 2 ) K af3{Qo, q)- 

(12) 

This expression is valid for arbitrary spectrum of magnetic turbulence including anisotropic 
distributions. 



-6 - 



We consider here the DSR produced by accelerated nonrelativistic electrons interacting 
with the MHD turbulence. In MHD waves E ~ (v a /c)B, where v a is the Alfven speed, therefore 
the magnetic part of the Lorenz force is larger than the electric part for all electrons with v > v a . 
Assuming this condition to be fulfilled, we calculate only the DSR related to the magnetic field 
of the MHD turbulence; inclusion of electric field effect will further increase the DSR intensity. 

Since we are interested in overall spectral shapes and flux level of the DSR, rather than 
model-dependent details of the emission, we consider here the simplest case of the isotropic 
MHD turbulence: 



K aP = \ Uap - q -^A K(q)S(q - q (q)). (13) 

As we assumed v > v a , i.e., the electrons move faster than the waves, we can adopt the MHD 
turbulence to be quasi static, qo(q) = 0. When the MHD turbulence is isotropic, the accelerated 
electrons are isotropic as well, and so the radiation produced is also isotropic. Thus, we consider 
further the radiation produced into the full solid angle 



/njn=^- ? ( W ), (14) 
61XC 



where, like in the ultrarelativistic case jFleishmanll2006h . we introduce the scattering rate of 



the nonrelativistic particle by MHD turbulence q(u>): 



gM = lGI) 2 ^ / *(q)*(" + qv)dq. (is) 



To proceed further we have to specify the shape of the turbulence spectrum if(q); we 
adopt a single power-law down to the smallest (resonant to thermal electrons) scales: 

K ^ = ^ ^ = ^K- l (Bl), (16) 

where k = 27r/L with L the largest turbulence scale, (B^ t ) is the mean square of the turbulent 
magnetic field, and v is the turbulence spectral index. 

Then, substituting (fl6|) into ([15]) . integrating over cfq, 



/ dqK(q)5(uj + qv) = 2vr / dcos9- dq^6(u + qvcos9) = ^ J ^ = 



V V \ulJ V J J \Vpe V 
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where 



v pe = 6.74xlOV T e 



is the thermal velocity of the plasma electrons, Q(x) is the step function, and using the electron 
charge e and mass m for Q and M, we find 

?h = ^4^1 c-y fi - f^v v - - V (19) 

yV ' 2Z/ m 2 C 4 W V \UpevJ J \V p e Vj' V ' 

so the DSR spectrum produced by accelerated electrons reads 

8e 2 



e-q{u). (20) 



37TC 

Now we calculate the DSR power from N electrons with the spectrum described by Eq. 

oo 

P. = J LN(E)dE. (21) 

In fact, we are interested in the radio flux observed at the Earth. To transform this radiation 
power into the flux observed at the Earth, we change the variable u = 2irf, so that If = 
Then, the flux is 

, toW = IV? . 10 19 f (22) 

f 4ttRI u 2Rl u y } 

where R au = 1 au= 1.49 x 10 13 cm is the distance from the Earth to the Sun. 

To evaluate the DSR from the acceleration region of a solar flare, we adopt some typically 
assumed parameters of the acceleration site as follows: (a) the size of the site L ~ 10 8 cm; (b) 
the thermal electron number density n e ~ 10 10 cm -3 ; (c) the electron temperature T e ~ 10 6 K; 
(d) the energy density of the magnetic turbulence W st = ~ 10 3 erg/cm 3 . Accordingly, the 
total energy, W to t ~ W st L 3 ~ 10 27 ergs, corresponds to a very modest solar flare. We assume 
that the power-law tail of the accelerated electrons grows from Eq = (4 — 6)kT e with n(> Eq) 
specified by matching condition (J2J) and as the acceleration has started, the power law index 8 
changes from 8 ~ 3 while the break energy Eb r increases from 50 ~ 500 keV. 

Figure [2] presents the sequence of calculated DSR spectra for 11 different 8 values from 
7 ~ 3; the spectra are calculated for three different v values and for two different Eq values. 
The blue curves indicate larger 8, while the red ones show smaller 8. Then, Figure [S^a) presents 
the DSR spectra for three different temperature values, T = (1, 3, 10) • 10 6 i^. In addition to 
spectrum shapes, light curves of the radiation at different frequencies can be informative. To 
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Fig. 2. — The DSR total flux density spectra calculated for the observations from the Earth 
surface for Eq = 4kT (upper panel) and Eq = 6kT (lower panel) for 11 different 5 from 7 ~ 3 
and three different v, indicated in the panels. The blue curves indicate larger 5, while the red 
ones show smaller 5. 



estimate the light curve behavior we adopt a soft-hard-soft sp ectrum evolution as follows from 
theory of spectrum evolution for the stochastic acceleration (IBykov fc Fleishman! 120091 ) . and 
which is typical for impulsive flares, with the electron energy spectral index 5(t) changing from 
8 to 3 and then back to 8, while the break energy E^ increasing all the way from 50 keV to 
~ 1 MeV. Figure El^b) shows the corresponding model light curves at a few frequencies around 
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(a) (b) 

Fig. 3. — (a) DSR spectra for three sources with 5 = 7, E = AkT e , and different electron 
temperatures, T e = (1, 3, 10) • 10 6 K; (b) DSR light curves for the case E = AkT ei v = 1.2, 
E br = (E br (t = 0) + ^ • t) keV, where E^t = 0) = 50 keV, ^g* = 45 keV/s, and soft-hard- 
soft evolution of the accelerated electron spectrum, 5 = 5 max + a ■ (t — t ) 2 , where <5 max = 3, 
to = 5 s, and a = 0.2 s~ 2 . The apparent symmetry of the model light curves is provided by 
assumed symmetry of the spectral index variation and adopted no variation of MHD spectrum; 
relaxing any of these assumption will result in asymmetric light curves. All the light curves 
peak almost simultaneously, so no "cross-correlation" delay is expected, although onset of 
higher frequency light curves is delayed, thus, the total duration of higher frequency emission 
is shorter. 



the spectrum peak. One can note from the figure that higher frequency light curves have a 
somewhat shorter duration, although peaking at the same time; so no appreciable time delay 
between the light curves is expected. 

We note that the DSR spectra are very narrowband, much narrower that typical gyrosyn- 
chrotron spectra. The high frequency slope of the DSR spectrum can easily be evaluated from 
Eqns. , (EH), (USD, and (EDI), F f oc f~ 2S . Thus, the DSR high frequency spectral index 
varies from 11 to 3 as the spectral index of accelerated electrons changes from 7 to 3, while the 
GS spectral index would vary from 5 to 1 for the same range of 5 variation. The peak flux of 
the DSR is highly sensitive to the turbulence spectral index (specified eventually by the MHD 
cascading law), while less sensitive to the plasma temperature and electron spectral index. 
The peak flux can be very large (up to a few hundred sfu), which makes it easily observable 
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even by full sun radio instruments. If so, the corresponding radio emission must have been 
widely observed by available radio spectrometers working in the decimetric and/or microwave 
range. Indeed, there is a class of radio bursts with the properties resembling the DSR proper- 
ties described here — it is the class of narrowband decimetric and microwave continuum bursts 
(including type IVdm), which, we suggest, may contain burst-candidates to the radio emission 
from the regions of stochastic acceleration in solar flares. Although this interpretation is tempt- 
ing, spatially resolved radio observations of the DSR will be needed to confirm it, to locate the 
region of stochastic acceleration, and study it in detail. Another plausible candidate for radio 
emission from stochastic accel eration episodes i s so-called transient brightenings, whose radio 
spectra are often narrowband (IGary et al.l 119971 ). 



3. Gyrosynchrotron Radio Emission from a Collapsing Trap 

Let us consider another model, a collapsing magnetic trap, which can efficiently accelerate 
charged particles. Unlike the stochastic acceleration models, no turbulence spectrum is essential 
to accelerate particles in the collapsing trap model; however, some spectrum of 'pre-accelerated' 
particles is needed, otherwise, the collapse of the trap will only give rise to plasma heating 
without nonthermal particle generation. 

Therefore, we assume that just before collapsing the trap contained both thermal plasma 
and nonthermal electron population with a power-law spectrum. To be specific, we adopt the 
initial conditions as follows: (a) the magnetic field strength Bq = 30 G; (b) the minimum and 
the maximum energy of the power-law spectrum E min = 0.01 MeV, E max = 1 MeV; (c) the 
thermal electron density n t h = 10 9 cm -3 and the non-thermal electron density n r i = 10 7 cm -3 ; 
(d) the source size L = 10" 



During the trap contraction, t he number of acceler a ted el ectrons evolves. For our modeling 



we adopt a solution obtained by iBogachev fc Somovl (120051 ). see Figure HI which takes into 



account the betatron and Fermi acceleration and the particle escape from the trap via the loss 
cone: 



where 



N = N - lVK ^ b (23) 

y/1 + (b m - 1)P 



b = b{t) = B(t)/B (24) 
l = l{t) = L(t)/L (25) 
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so b(t) changes from 6(0) = 1 to b m — B m /B , B m is the largest magnetic field value at the 
end of the trap collapse, and l(t) deceases from 1(0) = 1 to a very low value, say, l(t max ) = 0.1. 
For the sake of simplicity we assume a self-similar contraction of the collapsing magnetic trap. 
In this case, evolution of all parameters of the trap is uniquely defined by their initial values 
and the dimensionless source scale l(t). Thus, for any given contraction law, l(t), we can 
easily calculate the corresponding time history of all other relevant source parameters, as the 
magnetic field, the thermal electron number density, the source volume and the projected area, 
and the evolution of the nonthermal electron spectrum (jBogachev fc Somovl 120051 . 120071 ). For 
our modeling we adopt that the trap volume V linearly decreases with time during the trap 
contraction from 10" 3 to l" 3 ; we adopt 10 s for the trap collapse time, which is a few Alfven 
times (r a ~ L/v a ) for the trap parameters used. 




(a) (b) 

Fig. 4. — Left: The number of electrons in the trap according to the lBogachev fc Somovl (120051 ) 
model. Right: GS spectra in frequency range 0.1 ~ 10 4 GHz; the solid spectra are plotted for 
t — 0, 1, 9 s; the dashed-doted lines correspond to the end of collapse at t — 9.75, 9.975, 
9.9975 s. In fact, we take into account both GS and free-free contributions to the radiation 
intensity, however, the free-free component is inessential during the loop collapse: it is only 
noticeable at high frequency flat regions of the early spectra. 

Thus, we can straightforwardly calculate the GS spectra at different time frames and the 
radio light curves at different frequencies within the adopted collapsing trap model. Figure H] 
displays the GS spectra at different moments of the trap contraction. In agreement with a 
statement made in the previous section, at initial phase of acceleration the GS flux is small 
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(less than 1 sfu), which can only be recorded by high sensitivity spatially resolved observations. 
However, during the trap contraction the magnetic field increases and the fast electrons are 
accelerated, which all together lead to a significant increase of the peak flux and the peak 
frequency of the radio emission produced at the acceleration site; thus the radio emission 
becomes easily detectable by available radio instrument soon after the trap starts to contract. 




(a) 



(b) 



Fig. 5. — The light curves (a) and the normalized light curves (b) of the GS emission from a 
collapsing trap at a number of fixed frequencies 5, 10, 17, 34 GHz. The initial trap parameters 
and the contraction law are described in the text. A significant time delay between the light 
curves, comparable to duration of the light curves, is evident from the Figure. 

Figure [5] presents the light curves of the emission at a number of fixed frequencies, 5, 10, 
17, 34 GHz. Within the adopted model the peak flux increases with frequency, see Figures HJ [5] 
in fact, this increase ma y become less pronounced if the coulomb losses in the collapsing trap 
are taken into account (IBogachev fc Somovl 120091 ). A distinctive feature of the light curves, 
contrasting to that of DSR produced from the stochastic acceleration sites, is a noticeable time 
delay: the higher frequency light curves are delayed rela tive to lower frequency lig ht curves; this 
time delay will be present even when the coulomb losses (IBogachev fc Somovll2009l ) are included. 
A time delay in the sense predicted by our modelin g is frequently observe d in solar flares, in 
particular, in those with quasiperiodic pulsations ([Fleishman et al.l 120081 ) . Observationally, 
however, the GS emission from a collapsing trap can be contaminated by GS emission from 
trapped electrons produced by previous acceleration episodes, so unambiguous detection of the 
GS emission from a collapsing trap itself requires additional accurate analysis to separate the 
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contributions, which as yet has not been performed. 



Discussion 



There are many models in which electrons can be accelerated to nonthermal energies. Some 
mechanisms accelerate a tiny fraction of the electrons, which can only be observed via coherent 
radio emissions (e.g., type III bursts produced by electron beams, or accompanying metric 
spikes), others produce more powerful acceleration, sufficient to generate observable incoherent 
radio emission from either the acceleration site itself of from a remote 'radiation site'. 

The idea of using radio obse rvations to probe en ergy release/acceleration regions in flares 
has been around for awhile (e.g., iBastian et al. Il99sh. how ever, the studies focused mainly on 
coherent decimeter radio bursts. For example, iBena (119861 ) argued that decimeter narrowband 
millisecond radio spike clusters can be a signature of electron acceleration in flares, and, if so, 
the flare energy release must have been highly fragmented with e ach spike indicating a single 
energy release/acceleration episode. However, it has been found (lAschwanden &: Glidell Il992h 
that the radio spikes are frequently delayed compared with associated hard X-ray emission, 
implying the spikes ar e a secondary phenomenon associated wi th flares. Moreover, spatially 
resolved observations (IBenz et al.l 120021 ; iBattaglia &: Benzl |2009j) show that the spike sources 



are typically far away from main f 



radio spikes ([Fleishman et al 



main flare location (IGary fc Naqvi 



2003 



are locations. Even t hough higher frequency microwave 



Rozhansky et al.l 120081 ) can be produced at or around the 



20091 ). it seems doubtful that the coherent radio burst origi 



nate from elementary acceleration episodes (iFleishman fc Melnikovill998l ; iFleishman et al.ll2003 
Rozhansky et al.ll2008l : IBattaglia & Benzl l2009h . 



In contrast, in this Letter we have calculated incoherent radio emission from the accelera- 
tion region of a solar flare within two distinct acceleration models — stochastic acceleration by 
cascading MHD turbulence and regular (betatron and Fermi) acceleration in a collapsing trap. 
We have demonstrated that the radio emissions produced within these two competing acceler- 
ation models are distinctly different, which potentially allows distinguishing between them by 
the radio observations. In particular, we have found that the stochastic acceleration process is 
accompanied by a very narrowband DSR continuum radio emission, whose predicted properties 
are generally consistent with observed properties of narrowband microwave or decimetric (type 
IVdm) continuum bursts, thus, we suggest that some of those bursts can be produced from the 
sites of stochastic acceleration. 
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